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We present a method that permits the calculation of the dynamical correlation functions for 
r — ■ quantum systems. These are obtained by evaluating the generating functionals of the static moments 

Q\ ' of the relaxation functions in a self-consistent approximation that can be obtained by an appropriate 

^\ , extension of the effective potential theory. 
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The thermodynamics of quantum systems has been widely studied by the effective potential theory Jj]||]. This 

treatment has proved to be very accurate to determine equilibrium properties |3|-|fj], but encounters quite a number 

of theoretical and practical difficulties in being applied to dynamical phenomena. An extension of the formalism 

uniquely based on ingenuity does not allow for a serious control of the results and slips easily in ad hoc or even 

erroneous assumptions. On the other hand it is self-evident that an adequate way for calculating dynamical correlations 

I ■ would be highly welcome. It is our purpose to try to fill this gap and to provide an effective method to approach 

dynamical calculations in terms of path integrals. We shall mainly be concerned with the Kubo relaxation functions, 

"^ ■ that naturally appear in the framework of Mori theory [|5| . These functions are obtained by suitably defined scalar 

| t | products, the Mori products 

R A , B (t) = (A\B(t)) = du(A(0)B(t + iu)), (1) 
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where the observables A and B are taken such that (A) = (B) = 0, and braces denote the thermodynamic average. The 
Laplace transform of Eq. (Jlj), in turn, has a continued fraction representation whose coefficients are static normalized 
relaxation functions, i.e., time-independent normalized Mori products, of derivatives of the dynamical variables A(t) 
and B(t). In particular, the self-relaxation function So(4) = (FoIFq) -1 (Fo(t)\Fo) of an hermitian operator Fo(t) 

t-H , can be Laplace transformed and expanded in a continued fraction, namely Sj(z) = (z + 8j+i£j+i(z))~ , where 
6j+i = (FjlFj) -1 (Fj + i\Fj + i) and Fj denotes the so-called j-th fluctuating force. The quantities (Fj\Fj), can be 

iq ■ related to a combination of the first 2(j + 1) moments of the time series-expansion of Ho(t). While (Fj\Fj) with 
j y^ can be expressed in terms of static correlations, the quantity (Fq\Fo) requires the direct evaluation of the Mori 

C^ product. 

Experiments measure the spectral shape, given by the Fourier transform of the dynamic correlation function, related 
to E (z) by the "detailed balance" principle: 



1 



S(u) = {F \F ) —- -Sft(Z (z = iu;)). (2) 

g , e 

Therefore, the calculation ofS(u>) can be approached from the knowledge of the static quantities 5j up to a sufficiently 
large number j ' = J M , supported by some insight into the long-time behaviour of the continued fraction termination 

3 '(*)lli- 

We shall just provide an explicit formula for the direct evaluation of an arbitrary 5j , (Fq\Fq) included, by an effective 

potential approximation. We start from the generating functional in the hamiltonian path-integral form: 



Z[L,J] = j>V[x(u)} I V[p{u)] 
exp 



1 / / M 

— / du( — ip{u)x(u) + TC(p(u), x(u)) ~ TiL(u)p{u) — TiJ{u)x{u) j . (3) 

According to the effective potential method 0J3H we consider a quadratic trial action obtained by substituting 
Ti.(p(u),x(u)) with (l/2m) (p(u) — 1~]) 2 + (mui 2 /2) (x(u) — £) 2 + w, whose parameters w, m and uj depend on the average 
point (?7,£) = (/3h)~ 1 (fdup(u),fdux(u)) of each path. The effective Hamiltonian reads: 



H ct! ( V ,0=w(v,0+r 1 Hr 1 sinh/) 



(4) 



where / = f3hu>/2 rules the quantum character of the system. By defining the two-component vectors p = (?7, £) an d 
K{u) — t (L(u), J{u)), the approximated generating functional can be written as: 



Z [K] 
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du t P K(u) + - / 
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du 



/3h 



dv t K{u)$(u-v)K(v) 



(5) 



In Eq.@ we have introduced the 2x2 matrix $^ (u — v) with elements $n(u — v) = m 2 uj 2 & 22 (u — v) — to 2 cj 2 A/(u — v) 
and $12 (u — v) = — Q 2 i(u — v) = T/( u — v) , where 



Af(u-v) 



Tf(u — v) = im 



2mw sinh / 
dAf(u — v) 



cosh( \uj(u — v)\ — /) 



sinh / 



dv 



ih sinh( \u)(u — v)\ — f ) 
"2 sinh/ 



(u — v) — 6(v — u) 



(6) 



While Af is always well defined, the value of Tf for u = v is determined only if we specify the limit v — u — > ± 
that must be taken, reflecting the commutation relation of x and p at the same time. Moreover, a squared length 
a = h/(2mu>) (coth/ — 1//) and an energy E a — mui 2 a are naturally defined by the parameters of the system. It is 
also worth observing that A/(0) = Af((3h) = a and that A/(w — v) and T/(u — v) have a vanishing average in [0, (3K\. 
If we define the two-component vectors z = *(p, x) and y = (yi{u),y2(u)), after a lengthy but straightforward 
calculation, the following general formula can be derived in the low-coupling approximation [pi, i.e. lo independent of 
r\ and £: 
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(7) 



Here, TV = Z [0] is the normalizing factor, and the double braces denote the average over the Gaussian distribution 
of the variables yi v (u u ), defined by the second moments 



«w(f)yj(«)» = *«(«-«). 



(8) 



The expression (|7|) with the definition (|8|) represents the main result of our paper. Complicated static Mori products, 
i.e. moments of any order, can be evaluated by this last equation. Static correlations are also obtained performing 
the appropriate limit u — v — » 0. 

In order to proceed on our path, we first look for an explicit expression for averages of two functions separately 
depending on momentum or displacement. A direct calculation leads from Eq. (0) to the following expansion: 



T 2 u 



A {Zi (v))B {zj (u))])=Mj^e-WM 



exp 
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dydz 



A(y)B(z) 



$ij(u-v) 
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( z ~ Pi)} } 



(9) 



with &kk = ^fcfc(O). We note that, when static quantities are approached, i.e. the two imaginary times become equal, 
the particular ordering of operators depending on both p and x has to be specified. Any of the two possible orders, 
p-x or x-p, is obtained from the equation above by simply taking the appropriate limit v — u — > . 

The final step involved in the calculation of higher order moments is the integral expression for the average of the 
imaginary time ordered product of two operators A and B at two different times u and v, both depending on momenta 
and displacements. Without loss in generality we can consider p-left ordered operators, so that A — A\ (p(u + ))A2 (x(u)) 
and the similar for B. By specializing Eq.(ll) to this situation and with evident meaning of the notation, we obtain 



T u A(p(u),x(u))B(p(v),x(v)) ) = (T u i 1 (p( W +))i 2 ( a ;( M ))B 1 (p( U + ))B 2 ( a :(t;)) ) = 



= N f^ie-PKtMO 



2nh 



4 A 1 (y 1 )A 2 (y 2 ) B x (y 3 ) B 2 (y 4 ) 
& y 



(2 7 r) 2 det 1/2 (* l 



exp 



-9 2(W-/°*)*« 1 (» , i-Pj) 



(10) 



where now p — {t),£,7),£) is a four-component vector and we have introduced the 4x4 matrix W with elements 



vEr n = m I uj' I '$>22 = *33 = m 2 uj 2 ^ iA == m 2 uj 2 a, *i 2 = *2i : = *34 = *43 
m 2 w 2 *42 = m 2 u> 2 A f (u - v) and * i4 = -*23 = -*32 = *4i = T/(u - v). 
From Eq.(|lO|), when (u — v) — > + , we find the static average: 



-ih/2, *i 3 = m 2 W 2 *24 = ^31 



(AB) = N 



drjdl; 
2-nTi 



({ABfie-PH-vfa® 



(11) 



indeed recovering our previous result [§,gU|] . 

As far as the Mori product is concerned, the well known series expansion (2ne)^ 1 ^ 2 exp{— x 2 /(2e)} = 
zC^°=o (l/ n ( e /2) n (r 2 ™ - * 1 (x) appears to be an efficient tool to approximate the static Mori product of general op- 
erators, as well as their dynamical correlations, when the scales of the quantum fluctuations in the system, ruled by 
h and the natural length scale a, are small. Indeed, as the averages of A/(u) and Tf(u) in [0, /?7i] vanish, we have 



(A(p,x)\B(p 7 x))=N(Ph) y^ e -/^.«(r?,£){ ({A(t,,Q)) {{B(r,,£))) 

AqqB qq )<Tf + (A pp B qq — 



2m 2 ui 2 A pq B pq 



^^pq&pq 



+ AqqB pp )^f\ +o(a 2 ,?i 2 )}, (12) 



- (m 2 ^A pp B v 

where the subscripts of A and B denote the derivatives, while 07 and /i/ are the averages of A^ and F 2 * in [0, (3K\. 

Therefore, as <7/ is of order a 2 and \if of order h 2 , at the lowest order (A(p, x)\B(p, x)) reduces to the "classical like" 
average of the product of the Gaussian spreads of the two operators taken at the same order: 



(A(p,x)\B(p,x)) =Nph 



di] d^ 
2nh 



((A( v , 0)) ((5(77, 0» e _/3W -« ^ ® + o(a, h) . 



(13) 



In the following discussion, for the sake of simplicity, we consider a particle in a potential with a standard Hamil- 
tonian 



H(p, x) = p 2 /2m + V{x) ; V(x) - V"(x ){x - x ) 2 /2 = gv(x). 



(14) 



where xo is the absolute minimum of the potential, and the quantum coupling constant g has been defined. The 
effective Hamiltonian is therefore Ti. a (?7, £) = rj 2 /(2m) + V fl (t;), where m is now the particle mass, i.e. a constant, 
and V^CO = ((V(0>) + 0' 1 ln(/- ie sinh/) . 

The simplest physically relevant quantity that can be discussed is given by the imaginary time ordered correlation 
of two dynamical variables A(x) and -B(x), dependent upon the coordinates only. For standard Hamiltonians the 
integration over momentum can be done and we get the following explicit expression for the correlation function: 



(T u A(x(v))B(x(u))\) =N 
2A{y)B(z 



2-k(3Ti 2 
expj 



(y - z) 2 



dy / dz 



A[a - Af(u - v)] 



J expj 



{y + z-2Q 2 ■ 
4[a + Af(u — v)]. 



47r(a - A f (u - v)) 
and 



1/2 



1/2 



(15) 



From the last equation, or by specializing equations (fl 
obtained: 



47r(a + Af(u - v)) 
the expansion for the Mori product is finally 



(i(x)|B(x)) = Nf3h 



2wf3h 2 



««e e -^-r(0|^)B(0 + I [A"(0B(0 + A(0B"(0 



A {lv) (0-8(0 + +A(£)£ M (0 



^A"(0B"(0]}+o(a 2 ). (16) 



Notice that for the displacement-displacement correlation function, as probed by neutron scattering, only the first 
term in the last equation (or in Eq. (|12j)) survives, and gives the total contribution to (j4(x)|-B(x)) in the effective 
potential approach. Comparing with Eq.(|ll|), we have (x 2 ) — (/3h)~ 1 (x\x) + a . This shows the twofold contribution 
of the quantum fluctuations, i.e. the spreads of the potential and of the quantity under averaging. The Mori product 
takes into account the spread of the potential only. A result similar to our Eq.(|l2|) was proposed in |10| with ad 
hoc assumptions and using only the first cumulant. As we can calculate the zero moment, (a;)*- ) = (x|x), and since 



the higher moments are given by the other previous equations, the dynamic correlations (x(t)x) can be eventually 
determined by the continued fraction expansion [mq] . 

We like to conclude with some remarks on another possible dynamical extension of our theory. We can assume that 
our system evolves with the effective Hamiltonian as found for the thermodynamic behaviour. This means that we take 
the commutator for its quadratic part, while we consider the Poisson brackets for the one-loop-renormalized nonlinear 
part. In order to do this, we have to use a representation of the operators describing the dynamical variables which 
permits to unify quantum and classical dynamics, as, for instance, the Weyl-representation |3lp|] or the Liouvillian 
scheme M. To better understand such point we recall in the first place that harmonic oscillators evolve by the same 
law both in classical and in quantum dynamics: the differences between quantum and classical statistical evolution 
are due to the thermal occupation numbers, that are static quantities. For real systems, the quantum deviations from 
the harmonic behaviour are ruled by the coupling constant <?, related to quantum non linearity. If g is vanishing, the 
non commutativity of operators at different times can be neglected and Eq. (fi~3) can be assumed to hold its validity at 
different real times, provided also that a is small enough, so that the Gaussian spreads of A and B do not overlap. At 
this level, the averages in time of the quantities ({A)) and {(B(t))), evolving with the effective potential V eB (£) |10|,|ff|, 
can provide an approximation for the time-dependent Mori product (A\B(t)). Finally its Fourier transform gives then 
the relaxation function Ra,b{u), that is connected to the Fourier transform Ca,b (w) of the dynamical correlation 
function (AB(t)) by means of the fluctuation-dissipation theorem , as in Eq.(2). This last operation simply restores 
the correct static quantum occupation numbers. 

We stress that the results become exact when g — > 0, i.e. for quantum harmonic oscillators with a classical nonlinear 
interaction term. Of course, the other exact limit is the classical system. For finite values of g, the validity of this 
scheme involves also the amplitudes of the Gaussian fluctuations ruled by the parameter a. Therefore, there is the 
same behaviour found for approaching static correlators with the effective Hamiltonian [p|Jl2| . In that case, the lowest 
limit of temperature was found to be related to both parameters g and a. For the dynamic correlators, this procedure 
yields a good approximation for times up to the order of h/3, for which the use of the effective potential makes sense 
in the calculation of the static quantities at lowest order, reproducing for instance a correct second moment for the 
displacement-displacement dynamic correlator with a well-behaved classical long time decay. 
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